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Abstract 

In this article we will present a study of the well-known Kardar-Parisi-Zhang(KPZ) model. Under 
certain conditions we have found analytic self-similar solutions for the underlying equation. The 
results are strongly related to the error functions. One and two spatial dimensions are considered 
with different kind of self-similar Ansatze. 

PACS numbers: 64.60.Ht, 02.30. Jr, 68.35.Fx 



1 



I. INTRODUCTION 



Growth of patterns in clusters and sodification fronts are a challenging problem from 
a long time. Basic knowledge of the roughness of growing crystalline facets has obvious 
technical applications The simplest nonlinear generalization of the ubiquitous diffusion 
equation is the the so called Kardar-Parisi-Zhang [2] model obtained from Langevin equation 

^ = vV 2 h+±(Vhf + r,(x,t) (1) 

where h stands for the level of the local profile. The first term on the right hand side describes 
relaxation of the interface be a surface tension which prefers a smooth surface. The second 
term is the lowest-order nonlinear term that can appear in the surface growth equation 
justified with the Eden model and originates from the tendency of the surface to locally grow 
normal to itself and is non-equilibrium in origin. The last term is a Langevin noise to mimic 
the stochastic nature of any growth process and has a Gaussian distribution usually. The 
right hand side eventually may contain a constant term corresponding to a driving force that 
generates an increase of the surface and which will be discussed later. There are numerous 
studies available about the KPZ equation in the literature in the last fifteen years. Without 
completeness we mention some of them. Hwa and Frey [3] investigated the KPZ model 
with the help of the self-consistent mode-coupling method and with renormalisation group- 
theory which is an exhaustive and sophisticated method which uses Green's functions. They 
considered various dynamical scaling form of C(x,t) = x~ 2x C(bx,b z t) for the correlation 
function. Lassig shows how the KPZ model can be derived and investigated with field 
theoretical approach 4]. In a topical review paper Kriecherbauer and Krug |J derives the 
KPZ model from hydrodynamical conservation equations with a general current density 
relation. 

One may find similar models, which may lead to similar equations to the one presented 



above, modeling the interface growth of bacterial colonies 



6]. There is a more general 



interface growing model based on the so-called Kuramoto-Sivashinsky equation which is 
basically the KPZ model with and extra — V 4 /i term on the right hand side of (1). In the 
following we will not consider this model, because the self-similar Ansatz cannot be applied 
to this equation. 
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II. METHOD FOR SOLUTION 



To obtain analytic solutions of the KPZ equation in the following we will use the self- 
similar Ansatz which can be found in 10] 



T{x,t) = t-°f :=t- a f(co 



(2) 



where T(x, t) can be an arbitrary variable of a PDE and t means time and x means spatial 
dependence. The similarity exponents a and /3 are of primary physical importance since 
a represents the rate of decay of the magnitude T(x,t), while (3 is the rate of spread (or 
contraction if /3 < ) of the space distribution as time goes on. The most powerful result of 
this Ansatz is the fundamental or Gaussian solution of the Fourier heat conduction equation 
(or for Fick's diffusion equation) with a = (3 = 1/2. (Exponents with a value of one half 
in another language means a typical random walk process.) These solutions are visualized 
on figure 1. for time-points t\ < £2- Solutions with integer exponents are called self-similar 
solutions of the first kind, non-integer exponents mean self-similar solutions of the second 
kind. The existence of self-similar solutions exclude the existence any kind of characteristic 
time scale as well. 



Applicability of this Ansatz is quite wide and comes up in various mechanical systems [SHlO] , 
transport phenomena like heat conduction rjj or even the three dimensional Navier-Stokes 



equation 



12|. 



In the following we will consider one spatial dependence of the KPZ equation. Calculating 
the first time and space derivative of the Ansatz ([2]) and pluging back to ([I]) (first we consider 
no noise term i](x,t) = 0) we get the following constrains for the exponents: a = and 
— 1/2. The remaining non-linear ordinary differential equation(ODE) reads 



u ^ PI/ \ 
2 + 2 nu) 



0. 



(3) 



The result is proportional with the logarithm of the error function 

2vln 



\c\yfwD er / '(2 V / P)]+C2 
2v 



A 



(4) 



where erf is the error function 14j and C\ and C2 are integration constants. The role of C2 is 
just a shift of the solution. From physical reasons v the surface tension should be larger than 
zero. From analysis of the solution Eq. fl4]) the value of A should be positive as well. Figure 
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FIG. 1: A self-similar solution of Eq. ([2]) for t\ < t<i. The presented curves are Gaussians for 
regular heat conduction. 
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FIG. 2: The self similar solution Eq. ([4]) for c\ = c% = 1. Solid line is for A = z/ = 1, dashed line 
is for v = 1 and A = 2 and the dotted line represents v = 2 and A = 1. 
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2. presents two solutions with c\ = ci = 1 and for three different A and v combinations. 
Note, that all solution has the same simple qualitative behavior, a ramp-up and a converged 
plateau. Figure 3. shows a three dimensional x,t dependent solution of the original equation 
for ci = C2 = A = z/ = l The function has a similar structure like Fig. 2. a ramp- up and a 
convergent plateau. 

Let's consider some analytic noise term for the 1+1 dimensional case. From the self- 
similar Ansatz (2) it is obvious that the noise term rj(x, t) should be a function of u = x/t^, 
on the other side rj(x,t) should be a distribution function as well. Therefore we first tried 
the following noise terms: exp(-u) and exp(—u 2 ). Unfortunately, no analytic results could 
be found in any closed form. For the 1/(1 + u> 2 ) (Lorenzian) noise term we got analytic 
results which can be expresses in sophisticated terms of Heun C special functions which we 
not mention here. A second analytic result was found for the 1/u 2 noise term which can 
be evaluated with a tedious combination of the Whitakker M and Whitakker W functions 
which we also skip here. 

If the surface growth is fed by some mechanism with a driving force, like the constant v 
term that may appear in addition on the r.h.s. of Eq. 1. and which may be suggested by 
the work of [l3|, then other analytical solutions become available which can be expressed 
via Kummer M and Kummer U functions. If the constant driving term in its magnitude is 
not v (not equal with the parameter of the Laplacean term) but a general e, the solutions 
are even more complicated. There is no general closed form for this case where all the three 
parameters are free. If we fix A = v = 1 than three independent cases are available. For 
< e < 1 the solution is like in Eq. 4. If e = 2 the solutions contains a complex part, 
however for large (e = 10) the solution is evaluated via the combination of hypergeometric 
functions. Note, that all these solutions contain two integrations constants c%, C2 which have 
non-linear dependence in the solutions. 

At last we investigate the two spatial dependent cases where we generalize the self-similar 
Ansatz. In our former studies it came out clearly that the self-similar Ansatz of (2) can be 



generalised in many ways 111] / 2: 



h(x,y,t)=r a f(^^ :=*"»/("), (5) 

where F[x, y] can be understood as a parametrization of a 2 dimensional curve. For F[x, y] we 
considered various functions like the most important \J 1 x 2 + y 2 = a which can be interpreted 
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FIG. 3: The self-similar solution of the original Equation for c\ = C'i = \ = v =\. 

as the usual distance or the L 2 Euclidean norm. Unfortunately, now for the recent KPZ 
model only the linear case = y — ax — b is available giving us the following ODE without 
any contradiction 



u(a 2 + !)/» + /V) 



| + ^V)(a 2 -l) 



0. 



(6) 



For the exponents the a = and (3 = 1/2 constrains are still valid. The general analytical 
solution is the following: 



X(a 2 - 1) 



In 



1 1 



A via 2 + r 



A 2 <^ Cx^atrf 



2 y v{a 2 + 1) 



v{a 2 + 1) 



a - c 2 




oj \ a — 



u(a 2 + l) . (7) 



Note, that this solution is very similar to the one dimensional one (4) there is a ramp-on 
and a plateau for < u> which means positive time and positive interface widths these 
are physical constraints of a real solution. Figure 4 presents a particular solution for the 
A = z/ = Ci = c 2 = l,a = 2 parameters. 
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FIG. 4: Solution of Eq. 7 for c x = c 2 = A = v = 1, a = 2. 

III. CONCLUSIONS 

In summary we can say, that with an appropriate change of variables applying the self- 
similar Ansatz one may obtain analytic solution for the KPZ equation for one or two dimen- 
sion sometimes even with some noise term. 
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